Dimer Mott Insulator in an Oxide Heterostructure 
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We study the problem of designing an artificial Mott insulator in a correlated oxide heterostruc- 
ture. We consider the extreme limit of quantum confinement based on ionic discontinuity doping, 
and argue that a unique dimer Mott insulator can be achieved for the case of a single SrO layer in a 
GdTiOs matrix. In the dimer Mott insulator, electrons are localized not to individual atoms but to 
bonding orbit als on molecular dimers formed across a bilayer of two Ti02 planes, and is analogous 
to the Mott insulating state of Hubbard ladders, studied in the 1990s. We verify the existence of the 
dimer Mott insulator through both ab initio and model Hamiltonian studies, and find for reasonable 
values of Hubbard U that it is stable and ferromagnetic with a clear bonding/anti-bonding splitting 
of order 0.65eV, and a significant smaller Mott gap whose size depends upon U. The combined 
effects of polar discontinuity, strong structural relaxation and electron correlations all contribute to 
the realization of this unique ground state. 



Recently, the growth techniques from semiconductor 
physics, such as Molecular Beam Epitaxy (MBE), have 
been increasing applied to transition metal and rare earth 
materials to create correlated heterostructurespJ The re- 
sulting atomic layer control promises the ability to de- 
sign orbital, spin, and charge states, create new emergent 
phenomena, and study fundamental physics of correlated 
quantum states in unprecedented new ways. A first step 
in this direction would be the creation of the simplest 
and most dramatic manifestion of electron-electron in- 
teractions: the formation of a Mott insulators, a system 
which would be a metal according to band theory, but in 
which instead electrons localize because their motion is 
jammed by their mutual short-range Coulomb repulsion. 

Mott insulators occur only when the electron density 
is commensurate with the underlying lattice, and typi- 
cally an (odd) integer number of electrons per atom is 
required. A charge density of one electron per atom is 
enormous, reaching of order n2d ~ 7 x 10 14 cm -2 for a 
typical perovskite structure even if these electrons are 
confined to a single atomic layer. This provides a chal- 
lenge for heterostructures, as this n2d is already an order 
of magnitude larger than can be achieved in the high- 
est density semiconductors, and even if it is created, 
the electron density per atom will be greatly reduced by 
the electrons tendency to spread out in the third dimen- 
sion. In this paper, we show that these difficulties can 
be overcome by judicious design of a Dimer Mott Insula- 
tor (DMI), a state envisioned decades ago in the context 
of one-dimensional Hubbard ladders,^ and created here 
in two dimensions at a single monolayer of SrO embed- 
ded in a GdTiOs matrix. In the DMI, the requisite high 
charge density is achieved by combining ionic disconti- 
nuity doping, quantum confinement, and the formation 
of electronic dimers. The dimers are bonding orbitals on 
electron pairs, to which electrons are Mott localized in- 
stead of to individual atoms. Dimer formation halves the 
charge density needed to reach the Mott state, relative 



to the usual single atom localization, and is crucial to the 
success of our scheme. We combine ab initio and model 
calculations to establish the existence and nature of the 
DMI state theoretically. 

The starting point for our work is the polar/ionic dis- 
continuity, which induces a large net charge, typically 
half an electron per planar unit cell, at a polar to non- 
polar interface. The polar discontinuity has been iden- 
tified as a possible mechanism of doping in many oxide 
interface studies,^ but has only recently been quanti- 
tatively verified systematically 4 . In MBE grown het- 
erostructures of GdTiC>3 and SrTi03, a carrier density 
of ri2d = 3.5 x 10 14 cm -2 (=1/2 e~ per planar unit cell) 
for each GdTiOs /SrTi03 interface has been systemati- 
cally observed by Hall coefficient measurements 4 . These 
electrons fall into the empty d states of the SrTiOa, and 
consequently high 3d carrier density can be achieved by 
confinement in narrow quantum wells of SrTiOa embed- 
ded in thicker GdTi03, with n3^=n2^/w, where w is the 
well width. Recent transport experiments^ showed that 
such wells with a width of a few SrTi03 unit cells are 
indeed strongly correlated metals with ultra-high carrier 
density n^d = 7 x 10 14 cm -2 arising from two interfaces, 
corresponding to 1/2 + 1/2 = le~ per planar unit cell. 

To maximize the 3d electron density and approach 
the Mott limit, we take this approach to its logical end 
and consider the case of a single SrO layer embedded in 
GdTiOs. However, even in this case of ultimate confine- 
ment, we do not achieve a 3d density of 1 e~ per atom. 
This is because the doped electrons go symmetrically into 
the two interfaces Ti02 layers on either side of the SrO 
plane. While this situation appears unfavorable for the 
formation of a Mott insulator, all is not lost. On funda- 
mental grounds, the condition for the formation of a Mott 
state depends only on the charge per unit cell defined by 
the translational symmetry of the system. Here the two 
interface planes form a bilayer, with translational sym- 
metry only within the plane, and there is indeed a unit 
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FIG. 1: (Color online) One unit cell of the relaxed structure 
of the superlattice (SrTi03)i(GdTiOs)3. Ti06 octahedra are 
drawn in blue, to emphasize tilts.. The superlattice repeats 
periodically along the z direction. The bilayer consists of Ti02 
planes at z=0 and z=l. Here a and b denote the two different 
Ti sublattices. 



charge per planar unit cell. This indicates that a Mott 
state is possible in principle. To realize it, we must some- 
how induce the two Ti atoms in a unit cell to act as a 
single 'superatom'. 

It is instructive to view the bilayer on its side, with the 
Ti atoms projected into an x-z plane (here we use stan- 
dard cubic coordinates for the perovskite structure, and 
z is the growth direction). The Ti-O-Ti network then 
forms a ladder with Ti-O-Ti bonds between the two in- 
terfaces making the rungs of the ladder, and intra-plane 
bonds projecting to the ladders legs. In a Hubbard model 
description, symmetry dictates that the hopping ampli- 
tude for electrons along the legs, t, and along the rungs, 
are unequal. Correlated electron ladders were studied 
intensively in the 1990sPand in particular it was shown 
when t±/t is sufficiently large (approximately t±/t > 1 
for the one-dimensional Hubbard ladder), electrons form 
an unconventional Mott state of bonding orbitals on the 
rungs of the ladder - the one dimensional analog of the 
DMI. Qualitatively, we anticipate the same physics ap- 
plies to the Hubbard bilayers, provided t±/t is sufficiently 
large. 

Hopping parameters in transition metal oxides are 
largely controlled by the metal-oxygen-metal bond an- 
gle, due to the directionality of d and p orbitals, and 
are generally largest when the bond angle is closest to 
180°. Intuitively, we expect that the interlayer Ti-O- 
Ti bonds are the most SrTi03-like, while those within 
the Ti02 planes conform more closely to those of the 
GdTiOs. Since SrTi03 is nearly perfectly cubic, while 
GdTiOs is one of the most highly distorted titanates, 
this appears quite favorable. We checked this intuition 
with ab initio density functional theory (DFT) calcula- 
tions for periodic superlattices of the single SrO layer, 
(SrTi03)i(GdTi03) n , with n=3,5. Calculations were 
performed in the Wien2k 6 implementation and the gener- 
alized gradient approximation^ (GGA). An RKmax pa- 
rameter 7.0 was chosen with RMTs of 1.91 a.u., 1.69 
a.u., 2.29 a.u. and 2.27 a.u. for Ti, O, Sr and Gd, 
respectively. The DFT calculation is carried out in a 




4 



O z-coordinate 



FIG. 2: (Color online) Ti-O-Ti bond angles along (a) verti- 
cal direction and (b) in-plane direction with respect to the 
O (in the Ti-O-Ti bond) z positions. O along vertical direc- 
tion resides in the RO (R=Sr, Gd) plane, which is why the 
coordinate appears to be half-integer. 



y/2a x y/2a x c unit cell to allow for the possibility of 
octahedral tilts, where a is set to be the value of the ex- 
perimental SrTi03 lattice constant, 3.905A. The struc- 
tural optimization is done both on the atomic coordinates 
and c/a ratio, within the GGA + U approximation. We 
focus on U e ff=U-J=4eV on the Ti d orbitals, an accept- 
able value for the titanateiP. In addition, we further add 
U e ff=U-J=8.5eV on the Gd / orbitals since the energy of 
the occupied Gd / bands lie much lower than the Fermi 
energy in practice. This value will not affect the relevant 
electronic properties. 

We now discuss the key features of the relaxed struc- 
ture, which becomes independent of n for n > 3, and is 
shown in Fig. [T] First, relaxation is significant at the 
two interface layers, but decays quickly into the GdTiOs 
region. Second, the Ti-O-Ti bond angles are highly di- 
rection dependent near the interface, as shown in Fig. [2j 
The "vertical" Ti-O-Ti bond connecting the two inter- 
faces is slightly distorted, with a 160° angle. The next 
vertical Ti-O-Ti bond away from the interface is already 
highly distorted, with only a 3° angular difference from 
that of bulk GdTi03, and more distant bonds are nearly 
indistinguishable from bulk. On the other hand, the in- 
plane bonds are distorted but different from the bulk even 
in the interfacial layers, with bond angles of about 153°. 
Third, the Ti-0 bond length varies by only 6% between 
the longest and shortest bonds in the entire superlattice, 
and is less significant compared to the dramatic bond 
angle variations. 

Is the bond angle difference, together with other mi- 
nor structural relaxation effects sufficient to promote 
a DMI? We first check the electronic structure within 
the GGA+U approximation using the relaxed structure. 
Searching for possible magnetic structures, we obtained 
the lowest energy for a state with ferromagnetic align- 
ment of Ti spins within each Ti02 plane, with interfa- 
cial Ti spins antiparallel to those in the GdTiOs region 
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FIG. 3: (Color online) Layer-resolved electronic density of 
states of Ti 3d states of (a) sublattice a and of (b) sublattice 
b at z = 1 (at the interface); Ti of (c) sublattice a and of (d) 
sublattice b at z = 2 (center of GdTiOs region). The Fermi 
energy is set to zero. 



(energy difference compared to parallel to those in the 
GdTiOs region is small). With this configuration, and 
E/eff = 4eV, the interfacial and bulk density of states 
(DOS) is shown in Fig. [3j We observe a bulk gap of 
1.25eV, comparable to theoretical values in the literature, 
but remarkably a much reduced but still non-zero gap at 
the interface of approximately 0.2eV. This is a signature 
of the DMI state. Within GGA+U, the DMI persists for 
U e ff ^ 3.5eV, with both bulk and interfacial gaps reduced 
for smaller U . To see the DMI more directly, we decom- 
posed the DOS by Ti site and the 3 t<i g orbitals (denned 
by pseudo-cubic axes). Within the interface plane, the 
major DOS just below the Fermi level has no xy charac- 
ter, consisting instead of spin down predominantly xz/yz 
orbitals on alternating a/b sublattices. Consequently, 
we identify this state as the occupied bonding orbital 
of the DMI, with the xz/yz orbital degeneracy split by 
octahedral rotations. The anti-bonding state, centered 
at around 0.2eV, is again mainly composed of yz or xz 
orbitals. The separation (of subband centers) from the 
bonding state gives the bonding/anti-bonding splitting 
of 0.65eV. The gap is much smaller than this splitting, 
however, due the width of the bonding and anti-bonding 
bands. Away from the interface region, the electronic 
structure resembles that of bulk GdTiOs. For the Ti 
in the center of the GdTiOs region, we observe domi- 
nant xy/yz and xy/xz states alternating between two 
orthorhombic sublattices. Quantitative analysis shows 
the occupation at xy, yz and xz states are 37%, 16% and 
47% for Ti at sublattice a, similar to in bulk GdTiOs. 

To study the DMI state more explicitly, we constructed 
a model extended Hubbard Hamiltonian by extracting 
hopping parameters from the ab initio calculations. Us- 
ing the optimized structure obtained within GGA+U, a 
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TABLE I: Hopping parameters for two Ti02 interface layers 
from fits to the (SrTiOs)i (GdTiOs) 3 superlattice in unit of 
eV. The index, for example lay z, stands for sublattice a Ti 
at z = 1 with basis yz state. All the major hoppings (larger 
than 10% of the largest hopping magnitude) are kept here. 



pure GGA calculation was carried out. The eigenvalues 
of the / states of Gd were shifted manually away from 
the Fermi energy. We constructed 30 maximally local- 
ized Wannier function d 9 * 10 * (MLWF) around the Fermi 
energy, with the 3 ti g orbitals forming the basis. The 
hopping parameters were then calculated by evaluating 
the matrix element of the MLWF. Including these hop- 
ping parameters, the tight-binding model of the bilayer 
takes the form 

fftt = EE tT c L a c jna (i) 

(ij) mn,a 

where i and j are nearest neighbor sites, m and n are 
orbital indices and a is the spin index. The hopping 
parameters for the two interfacial Ti layers are tabulated 
in Table 1, which contains the full orbital dependence of 
the hopping terms. The maximum intra-plane and inter- 
plane hopping matrix elements are t ~ 0.35eV, t_L ~ 
0.63eV, respectively, which gives strong support for the 
DMI picture. The separation energy between bonding 
and anti-bonding state is consistent with this magnitude 
inter-plane hopping. 

We supplement these hopping parameter s w ith 
interactions to form the effective HamiltoniarP^ 

H = H tb + H int 

Hint = ^ ^ Minimi + U' ^ ^imt n in| 

(2) 

where U and U' represent on-site intra-orbital and inter- 
orbital Coulomb repulsion between up and down spin, 
respectively, and J is the Hund coupling. We will not 
restrict the condition to Slater-Kanamori interaction pa- 
rameter U' = U — 2 J and J = J' , but rather simply 
assume J = J' and explore the phase diagram by vary- 
ing all the other parameters. 
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FIG. 4: (Color online) Hartree-Fock phase diagram for Hund's 
couping J =0.6eV. Here, we abbreviate the phases as follows: 
PM+M = paramagnetic metal; CO+M= Weakly charge or- 
dered metal; CO+I = charge ordered insulator; FM+M = fer- 
romagnetic metal; FM+I = ferromagnetic insulator; AFM+I 
= antiferromagnetic insulator. 

Fig. [4] shows the phase diagram as a function of U and 
U' obtained by the Hartree-Fock approximation, with 
fixed Hund's coupling J =0.6eV. When the Coulomb re- 
pulsion is small, the ground state is just paramagnetic 
metal. In the unphysical region where inter-orbital re- 
pulsion is dominant, U' > U electrons reside in a single 
orbital per Ti and in fact charge order. Of most interest is 
the lower right part of the phase diagram, where U > U' . 
In this region, ferromagnetism arises for sufficiently large 
U. In the weak inter-orbital repulsion limit, the elec- 
trons are distributed in all three t2 g orbitals. Under this 
condition, the system is always metallic or semi-metallic, 
since those electrons in xy orbitals are non-bonding be- 
tween layers and hence metallic. For sufficiently large U' , 
however, the inter-orbital repulsion eventually disfavors 
and empties the xy states, in favor of the xz/yz orbitals 
which have lower energy through inter-layer hopping. In 
this way the bonding state becomes fully occupied and 
a (ferromagnetic) DMI is achieved. Further increase of 
U' enhances the gap and eventually prefers an antifer- 
romagnetic DMI, since Hund's coupling becomes ineffec- 
tive if there is strictly one electron per site, and super- 
exchange becomes dominant. For completeness, we also 
have tested J=0,leV, which shows the Hund's coupling 



enhances the ferromagnetic insulating state. For our best 
guess at physically appropriate values, e.g. U = 4.5eV 
and the Slater-Kanamori U' = U — 2 J, the ferromagnetic 
DMI obtains. The Hartree-Fock results for these values 
for the band gap and orbital ordering are very similar to 
those of the previously discussed GGA+U calculations. 

In summary, we have argued for the existence of a 
Dimer Mott Insulator (DMI) for a single SrO layer em- 
bedded in a thick GdTiOa matrix, using both ab ini- 
tio and model calculations. The DMI state is unique 
to the bilayer Ti02 structure created by a single SrO 
layers; we have indeed verified that a metallic state is 
obtained in GGA+U for the case of two SrO layers em- 
bedded in GdTiOs (see supplemental material). Insu- 
lating behavior in such structures experimentally should 
therefore be attributed to the combined effects of disor- 
der (e.g. SrTi03 thickness fluctuations) and interactions. 
The DMI for a single SrO layer could be experimentally 
probed by many experiments, including transport, op- 
tical measurements of the gap and bonding-antibonding 
splitting, and angle-resolved photoemission. Observing 
the magnetic structure is more difficult, but might be 
possible with ferromagnetic resonance or optical dichro- 
ism. This work suggests many directions for future the- 
oretical and experimental research. We anticipate that 
the Mott state can be controlled and modified by vary- 
ing composition, strain and the growth direction. It may 
be possible to create antiferromagnetic DMIs by varying 
the rare earth ion; however, theory is needed to gauge 
whether this also may destabilize the dimer formation it- 
self. Choice of substrate also effects the strain and growth 
direction of the titanate films. Whether the DMI persists 
when the GdTiOa grows along the (110) direction is an 
important question for future study. More speculatively, 
we might contemplate the possibility of superconductiv- 
ity induced by doping the DMI, by analogy to the super- 
conductivity predicted theoretically and observed exper- 
imentally in ladder systems. 
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